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ABSTRACT 



The parameter imbedding method for the Fredholm equation converts it 
into an initial value problem in its parameter A. We establish the method 
for general operator equations of the form [/ + f{X)]ip = 4>. It is particu- 
larly useful for studying spontaneous symmetry breaking problems, such as 
contained in the nonlinear Schwinger-Dyson equation. 
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I. Introduction 

A very useful technique in the solution of the Predholm integral equation 
is the parameter imbedding method.^ Consider the equation 

ilj{x) = xj K{x,y)i>{y)dy+ (j){x). (1) 

It is well-known that for = 0, the equation has nontrivial solution only for 
certain eigenvalues A = Aj, and that for A /=Aj, the solution is given by 

iP{x) = <P{x) + Xj R{x, y, X)(f>{y)dy, (2) 

where the resolvent function can be written in the form 

R{x,y,X)^D{x,y,X)/d{X). (3) 

The parameter imbedding method consists of a pair of differential equations 
relating D{x,y,X) and d{X), 

d\X) = 4^d{X) = - / D{x, X, X)dx (4) 



^D{x, y, X) - D{x, y,X)^ + -^ I D{x, z, X)D{z, y, X)dz, (5) 

d(0) = 1, (6) 



d{X) d{X) 

with the initial conditions 



D{x,y,Q) = K{x,y). (7) 

These equations can be integrated numerically with respect to A, yielding 
the solution to the original integral equation. Note that the eigenvalues of 
the original equation are just the zeroes of rf(A), also, the eigenfunctions are 
proportional to D{x,y*, X), for any fixed value of y = y*. The fact that 
d{X) and D{x,y,X) are analytic facilitate the numerical solutions greatly. 
Thus, the parameter imbedding method turns the original integral equation, 
with variable x and y, into an initial value problem with variable A. The 
latter is often much simpler and explicit solutions can be readily obtained. 
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For a practical application, the reader is referred to the solutions of the BS 
equation obtained earlier.^ 

However, for a number of problems in mathematical physics^, one encoun- 
ters integral equations which are not of the Fredholm type. For instance, 
a question of considerable current interest concerns systems which exhibit 
spontaneous symmetry breaking. They are often described by nonlinear, 
Hammerstein equations which are of the form 

^(x, A) = A 1 K{x, y)Fiy, ^(y, \))dy, (8) 

where F is a nonlinear function of tjj. The dependence on the parameter A is 
thus rather complex. In fact, it typically has bifurcation solutions which sig- 
nal the onset of spontaneous symmetry breaking. One may attempt to solve 
these equations by linearization and then apply the parameter imbedding 
method. However, because of the complicated A dependence, the imbedding 
method for Fredholm equation can not be used. 

The purpose of this paper is to find a generalization of the parameter 
imbedding method which can be used in this situation. To this end we have 
studied a general class of operator equations in a Hilbert space. They are of 
the form 

[/ + /(A)]V' = 0. (9) 

where we will assume that the operators /(A) and -^f{X) are linear, analytic 
in A, and of trace class"''^. (An operator A is of trace class if tr\A\ < oo). 
Provided we know the inverse operator, [/ -|- /(A)]~^, the solution would be 
iP=[I + /(A)]-V- Let us write 

R{X)^[I + f{X)r = ^^D{X), (10) 

where 

d{\) = det[I + fiX)]. (11) 

Then, it is found that the functions d{X) and -D(A) are analytic in A and they 
obey the following pair of parameter imbedding equations 

A,(A)=T,.|A;(A).i,(A)l. (12) 
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Further, we can determine the initial conditions for d{X) and D{X) at an 
arbitrary point A = Aq. Thus, one can integrate the imbedding equations 
and arrive at the solution R{X) for any A. The imbedding equations, Eqs. 
(12) and (13), are very simple, yet rather general. It should be applicable to 
a wide class of problems. 

We now comment on the background materials to be used in this work. 
Since we need to cite results in functional analysis extensively, it is difficult to 
make our paper self-contained. For this reason we will follow closely the four 
volume treatise "Methods of Mathematical Physics", by Reed and Simon^, 
for notations, concepts, and theorems. The only exception is that, in this 
paper, all operators will be denoted with a "hat", as in /(A), -D(A), etc. 

This paper is organized as follows. The analytic Fredholm theorem in 
functional analysis will be presented in Sec. II. In Sec. Ill and Sec. IV, 
we derive the two parameter imbedding equations, which are the central 
results of this work. The initial conditions accompanying these equations are 
determined in Sec. V, and concluding remarks are offered in Sec. VI. To make 
connections between our operator formalism and that of the usual Fredholm 
method, we have shown, in Appendix A, a specific realization which leads to 
the familiar results. Finally, we present in Appendix B an explicit formula 
for the partial trace, which also establishes its relation with the Plemelj and 
Smithies formula. 

II. The Fredholm Theorem 

In functional analysis, the Fredholm integral equation has been abstracted 
so that it is applicable to the general operator theory in a Hilbert space. We 
now quote the analytic Fredholm theorem which will be useful for this work.^ 
Theorem: Let D be an open connected subset of C. Let / = D — > C{TC) be 
an analytic operator-value function such that /(A) is compact for each XeD. 
Then, either 

A) [/ + /(A)]^^ exists for no XeD, or 

B) [7-|-/(A)]~^ exists for all X€D\S, where S is Bk discrete subset of D (i.e., a 
set which has no limit point in D). In this case, [/-|-/(A)]^^ is meromorphic 
in D, analytic in D\S, the residues at the poles are finite rank operators, 
and if XeS then —f{X)il) = ip has a nonzero solution in Ti. 

^From this theorem, we see that it is only necessary to ascertain that 
[/ + /(A)]~^ exists for some parameter value A = XgeD, then the operator 
/(A) would have the properties described in B). As we will show later, this 
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amounts to the condition that d{Xo) /=0 for some A = XoeD. Hereafter, we 
assume that this condition is always satisfied. 
Recall the definition, 

^(A) = [7 + /(A)]-^ (14) 

As long as XeD\S, we have 

^(A) + f{X)R{X) = /, (15) 

RiX) + RiX)fiX) = I. (16) 

Thus, consistency demands that /(A) and R{X) commute. Note that R{X) 
and /(A) are analytic in XeD\S and XeD, respectively. They and their 
derivatives are bounded in D\S. For later uses we also note that 

^^R{X) = -R{X)-^J{X)R{X), (17) 

which follows directly from Eqs.(15) and (16). 

III. The First Parameter Imbedding Equation 

We now prove our first parameter imbedding equation, Eq.(12). Recall 
that we demand both / (A) and ^/(A) to be of trace class. Since an operator 
which is of trace class must also be compact,^ the results of Sec. II are 
applicable to /(A) and ^/(A). 

Let's first define the determinant 

oo 

d{X) = detil + fix)) = Yl Tr[A\fiX))], (18) 

k=o 

where we used the definition^ 

Tr[A\A)]= (e^i ACi, A... Aei,,iei, A... Aie^J, (19) 

l<ii<...<oo 

Cj being a basis of the Hilbert space. Also, Tr[A''(74)] = 1. The series defining 
d{X) turns out to be uniformly convergent, which we now prove. 
Let us define the norm of a trace class operator, Aeli, by 

||i||i = ir|i|. (20) 
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This gives Xi a norm topology under which it is a Banach space satisfying 
the inequahty^ 

||A'=(i)|k<||i||tA!. (21) 

For any \eD,f{X) is of trace class, hence ||/||i is bounded. Let's write 
M = s'up||/(A)||i, for XeD. Then 

l|A'=(/(A))||i<||/(A)||JA!<MVfc! (22) 

It follows that the series for d{X), Eq.(18), is uniformly and absolutely con- 
vergent. We may thus differentiate term by term 

^diX) = Y.^Tr[A'^ifim- (23) 

k=o 

Or, using Eq.(19), and changing orders in the exterior products, we may 
write the result in the form 

^d{X) = J2k Y: (enA...Ae,„-/(A)e,A...A/(A)eJ (24) 

""^ k=o l<ii<...ik<oo 

From Appendix B, Eq.(B.l), this may be written as 

-d{X) = Y: kTr[-f{X) ■ Tr,_i(A^(/(A)))], (25) 

k=o 

where the partial trace operator is given by Eq.(B.6), 

k (_-]\m+l 

Tn.,[A\f{X)] = E ,^ \, r-^(A)rr [A'^-(/(A))] (26) 

m=l K -m+ 1 

Note that this equation for ^<i(A) is well-defined since ^/(A) is of trace class, 
rrfc_i[A''(/(A))] is bounded or of trace class, and the product of a bounded 
operator and one of trace class is also of trace class. ^ In addition, the series 
defined in Eq.(25) is absolutely and uniformly convergent, following a similar 
proof for the series of d{X), Eq.(18). We may thus interchange Tr and J2kJ 
resulting in 

d d ^ °° 

-d{X) = Tr[-f{X) E kTn.,[A'{f{X))]. (27) 

k=l 
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Introducing the definition 

D{X)^f:kTr,_^[A'{f{X))], (28) 

k=l 

we have 

^^(A) = Tr[^/(A)-^(A)]. (29) 

This is the first parameter imbedding equation. However, it remains to es- 
tabhsh the relation ^(A) = [/ + /(A)]"^ = D{X)/d{X), which we will do in 

the next section. 

IV. The Relation ^(A) = D{X)/d{X) and the Second Parameter 
Imbedding Equation. 

/,From Eqs. (15-16), in order to verify that R{X) — D{X)/d{X), it is 
sufficient to establish that 

D{X) + f{X)D{X)^d{X). (30) 



DiX) + D{X)fiX)=d{X). (31) 

Let us first note that D{X) and /(A) commute. This follows since the partial 
traces in D{X) are all polynomials in /(A). Thus, we need only to examine 
Eq. (30). Substituting the formulae for d{X) and D{X) from Eqs. (18) and 
(28), we have 

oo oo 

Y: kTn.,[A\f{X))] + f{X) Y: kTn.,[A\f{X))] 

fe=i ^ fc=i (32) 

= E Tr[A'-\f{X))]. 

k=l 

In this equation, we have used the convention Tr[A°(/(A))] = 1, /(A)° = 
1. Eq. (32) is vahd since, for any k, 

k Tn.,[A\f{X))] + Hm - i)Tn.,[A'^-\f{X))] 

= Tr[A^-H/(A))]. ^"'"'^ 

This follows from the Plcmelj-Smithies formula, Eq. (B.7), when we evaluate 
the determinant by expanding along the first row and using Eq. (B.9). Note 
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also that the special case, k=l, amounts to Tro{f{X)) — 1, which follows 
from Eq. (B.6). 

Having thus verified the relation R{X) = D{X)/d{X), we can now pro- 
ceed to derive the second parameter imbedding equation. This we obtain by 
substituting ^(A) = D{\)/d{\) in Eq. (17), derived in Sec. II. It follows 
immediately that 

which is the second imbedding equation. Further, we shall now deduce that 
D{\) is analytic for all XeD\S. It is already known that d{\) is analytic for 
XeD. Also, if d{X) ^0, the zeroes of d{X) are isolated. Prom the Predholm 
theorem in Sec. II, -R(A) = [/ + /(A)]~^ exists for XeD\S and is a mero- 
morphic function, whose residues at its poles are operators of finite rank. 
Since we already know that D{X) = R{X)d{X) is analytic in D\S, to prove 
its analyticity in all of D, we need only show that the rank of its zero of d{X) 
is greater than or equal to that of the pole of ^(A). For this purpose we may 
introduce the operator Px^ which projects Ti into the vector space spanned 
by the eigenvectors associated with X^eS. This vector space, {Ao}, is finite 
dimensional and has dimension dimlXo}. We decompose R{X) with respect 
to and (1- Pa„), 

^(A) = [/ + /(A) Aj-^ A„ + [/ + /(A)(l - Aj]-^(1 - A„). (35) 

At A = Ao, the second term is nonsingular, while for the first term, /(A)Pao 
is an operator in the finite dimensional vector space {Ao}. The inverse oper- 
ator [I + /(A)Pao]~^ can thus be obtained with the usual finite-dimensional 
techniques and is given by the ratio of a co-factor and the determinant of 
[/ + /(A)Pao]- The latter is a polynomial of degree no higher than dim {Ao}, 
so that the rank of its zero must be < dim {Ao}. Thus, we have proved that 
D{X) — R{X)d{X) is an analytic function in all of D (including S). 

V. The Initial Values of ri(A) and D{X). 

The parameter imbedding equations, given by Eq.(12) and (13), can be 
integrated provided we have the initial values of d{X) and -D(A) at some point, 
say A = Ao. This problem becomes trivial if it happens that /(Ao) = 0. In 
this case R{Xo) = [i" + f{Xo)]~^ = / so that d{Xo) = 1 and D{Xo) = /. In 
general, we can enforce this condition by introducing another parameter. 
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As long as d{X) ^0, we can always find a Ao so that d{Xo) /=0. Let us 
consider a new operator equation [/ + ^f{Xo)\^ = 4>, where ^ is a complex 
parameter. The operator-valued functions ^f{\o) and -^{^/{Xo)) = /(^o) 

are obviously analytic in ^ and of trace class. Defining [/ + C/(Ao)]~^ = 
Xo)/d{^, Ao), we have the imbedding equations with respect to ^: 

^d{^,Xo)^Tr[f{Xo)D{^,Xo)], (36) 

Ao) = Ao) - /(Ao)^(e, Ao)] (37) 

For these equations, the initial values at ^ = are obviously 

c?(0,Ao) = l, (38) 

^(0,Ao) = /. (39) 

We may therefore integrate the imbedding equations to the point ^ = 1, 
where 

rf(l,Ao) = rf(Ao), (40) 

D{l,Xo) = D{Xo). (41) 

Provided that d{Xo) /=0, these are then the initial values we need at A = Ao, 
which can be used with the original Eq.(12) and (13). 

VI. Concluding Remeirks 

In this work we have obtained a general method to solve the operator 
equation [I + f{X)]ip — 0. If /(A) and ^/(A) are hnear, analytic in A, and of 
trace class, then there is a pair of simple differential equations for the analytic 
functions D{X) and d{X), defined by ^(A) = [/ + /(A)]"^ = D{X)/d{X). The 
integration of these equations is rather straightforward so that one can obtain 
exphcitly the functions d{X) and D(X). 

Since the parameter imbedding method tackles a problem directly with 
respect to the parameter A, it is particularly useful when the solution ex- 
hibits intriguing properties in A. An example is the Hammerstein equation, 
which includes the familiar nonlinear Schwinger-Dyson (SD) equation. There 
are several advantages in the application of the imbedding method to this 
type of problems. First of all,the SD equation is known to have bifurcation 
solution. However, the criterion of bifurcation is contained in the function 
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d{X). As long as d{X) /=0, the implicit function theorem is valid and we 
have a unique solution. Thus, at each step of the integration of the parame- 
ter imbedding equations, the obtained value of d{X) enables us to determine 
whether there is a unique solution at A. Therefore, we can straightforwardly 
arrive at a solution for a wide range of the parameter A. Another nice feature 
of the imbedding method is that, at each step of the integration, the result 
obtained is directly the solution in its final form. It is thus very economical, 
computationally. Finally, the analyticity properties of d{\) and -D(A) are 
extremely useful, since we can always choose contours around possible singu- 
larities and integrate over smooth functions. All of the numerical calculations 
are therefore routine. 

These points are well illustrated in a concrete example, which we present 
in a separate paper where we have found numerical solutions to the SD 
equation. Our imbedding method enables one to obtain explicit bifurcation 
solutions to the SD equation. In the process we also uncovered some of its 
hitherto unknown solutions. 

The imbedding equations discussed in this work are very general. There 
seems to be little obstacle in applying them to a number of interesting prob- 
lems in mathematical physics. We hope to turn to these questions in the 
future. 

APPENDIX A: R(A) AND THE FREDHOLM RESOLVENT 

In this appendix, we examine the relation between the operator -R(A) 
and the usual Fredholm resolvent. They are not identical, which is why 
imbedding equations for the two cases also differ. The Fredholm equation 
(1) can be written as [/ -|- /(A)]V' = cf) with the definition 

/(A) = -A I K{x,y)^{y)dy. {Al) 

Its solution [Eq.(2)] may be written in the form 

= [/ + Ai?]0 (^2) 

with 

R(j>^ I R{x,y,X)(f>{y)dy. (A3) 
Thus, the relation between -R(A) = [/ + f{X)]~^ and R is 

I + XR^R{X)^[I + f{X)]-\ (A4) 
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or 

/(A) + Ai?[/ + /(A)] = 0. (A5) 

Thus 

it:=-^/(A)^(A). {A6) 

We may now proceed to find a relation of the operator imbedding equa- 
tions (12) and (13) for the Fredholm equations. Introducing the coordinate 
basis \x), so that {x\ip) — ip{x), and defining 

{x\f{X)\y)^-XK{x,y), (A7) 

{x\DiX)\y) = D{x,y,X), {A8) 

then Eq.(12) becomes 

^d(A) = Tr[^f{X)b{X)] =J J{x\ ^f{X) \y) {y\ D{X) \x) dxdy 

K{x,y)b{x,y,X)dxdy. {A9) 
Also, taking the matrix element of Eq.(13), we have 

d(X)—D {x, y, A) = D{x, y, X)—d{X)+ / D{x, x^ , X)K{x^ , z)D{z, y, X)dx^dz, 
oX dX J 

(AlO) 

but Eq.(A6) imphes that 

R{x,y,X) = j K{x,z)R{z,y,X)dz, {All) 

so that 

D{x, yA)^ j K{x, z)b{z, y, X)dz, {A12) 

Thus, Eq.(A9) is precisely Eq.(4). Also, Eq.(AlO) is equivalent to Eq.(5) by 
multiplying with K{x/ ^x) and integrating with respect to x. 

Although the two approaches are equivalent, the use of ^ = [/ + /(A)]~^ 
is convenient for general formulation. This is because the operator /(A) is 
compact, but in an infinite dimensional space, its inverse is not bounded. 
On the other hand, the identity operator is bounded, but not compact. The 
combination [7-|-/(A)]~^ strikes a balance and has nice invertibility properties 
similar to those of operators in infinite dimensional spaces. As for the Fred- 
holm resolvent, we note that from Eq.(A5), it satisfies I + XR[I + f~^] — 0, or 
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R = —[I + f ^] ^/X. The appearance of the unbounded operator / ^ makes 
R less useful compared to R. 

APPENDIX B: PARTIAL TRACE AND ITS EXPLICIT FORM 

Important properties of the partial trace will be examined in this ap- 
pendix. For this, we follow Simon, where the partial trace Trk-i is defined 
as an operator which projects operators in Ii{(S)''Ti.)to Ii{Ti). Its satisfies 
the relation 

Tr[CTrk-i{A\A))] = J] (e,, A ... A e^,, Ce,, A ie^, A ... A ie,J (51) 

l<ii<...ijt<oo 

for any C eIi{H). 

We now proceed to give an exphcit realization of the partial trace Trfc_i(A^(A)). 
Let us expand the right-hand side of Eq.(Bl) with respect to the first column 

^ (Cii A ... A ei^,Cei^ A ie^^ A ... A ie^J 

l<ii<...ifc<oo 

= - ^ A; [(eii,C'ejJ(ej2 A ... A ei^,iei2 A ... AicjJ 

ii=o l<?2<...<0O 

k-l 

+ IZ(-iyK>C'eiJ(ei,... Aei._, Ae^.^j A...,iei, A ...)] 

-1 oo _ 1 °o 

= T:E(en,Ce,jrr[A^-^(^)]-— ^ (e,^.,CeJ 

H=o ii=l l<i2---<oo 

{...ei._, A Ci.^,..., ...Aei._, A 1^.^,...). {B2) 

In deducing the first equality, note that the factor l/k\ in the definition 
of A^(yl). For the second equality, note that there is a factor (— l)''^-'^for 
bringing Aci^to the front of the exterior product. We may now expand the 
second term in Eq.(B2) into two terms, etc. Using the relations 

oo 

J2{e,„Ce,,)^Tr{C), (S3) 

oo oo 

E {e,.,Ce,,){ei„Ae,.) = {e,^ , C Aci.) = Tr{CA), (BA) 
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etc., we have 



Tr[CTr,^M\A)\^Y.i-^) 



m+l 



m=l 



{k — m)\ 
k — m 



It is now clear that the partial trace is given by the explicit formula 

k I -j^y((t+l 



\ k — m + 1 



A'^-^Tr [a"-"" (A)]. 



{B5) 



{B6) 



Note that this expression for the partial trace Tvk-i indeed projects jC{<Si'^H) 
into £(7i),as required. 

Finally, we show that this result for the partial trace is intimately related 
to the plemelj-Smithies formula, which is very important in the general 
Predholm theory. The plemelj-Smithies formula is obtained through the de- 
terminant 



A k-1 

Tr{A) k-2 



Tr{A''-^) ... 
Prom the following expression for the trace: 



k\ 



Tr{A) k-1 
Tr(i2) Tr{A) k-2 

Tr(i'') 






Tr{A) 







Tr(i) 



{B7) 



m 



we can write Eq. (B7) by expanding according to the first column 



m+l 



A;! 



m=l 



k — 171 + 



-A^'Tr = klATrk-iiA^A)]. {B9) 



Thus, the explicit formula of the partial trace (B6) is a generalization of the 
plemelj-Smithies formula. This result has been crucial in establishing the 
parameter imbedding theory for the trace class operator-valued functions. 
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